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Abstract 

The pressure, and the energy and entropy densities are determined for 
the SU (3) gauge theory in 2 + 1 dimensions from lattice Monte Carlo cal- 
culations in the interval 0.6 < T/T c < 15. The finite temperature lattices 
simulated have temporal extent N T = 2, 4, 6 and 8, and spatial volumes 
Ng such that the aspect ratio is Ns/N T = 8. To obtain the thermody- 
namical quantities, we calculate the averages of the temporal plaquettes 
P T and the spatial plaquettes P$ on these lattices. We also need the zero 
temperature averages of the plaquettes Po, calculated on symmetric lat- 
tices with N T = Ns- We discuss in detail the finite size (^-dependent) 
effects. These disappear exponentially. For the zero temperature lattices 
we find that the coefficient of Ns in the exponent is of the order of the 
glueball mass. On the finite temperature lattices it lies between the two 
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lowest screening masses. For the aspect ratio equal to eight, the system- 
atic errors coming from the finite size effects are much smaller than our 
statistical errors. We argue that in the continuum limit, at high enough 
temperature, the pressure can be parametrized by the very simple formula 
p = T 3 (a — bT c /T) where a and b are two constants. Using the thermo- 
dynamical identities for a large homogeneous system, this parametrization 
then determines the other thermodynamical variables in the same temper- 
ature range. 



1 Introduction 

The determination of the thermodynamical variables in the high temperature 
phase of QCD is one of the main goals of lattice gauge theory. The values of 
those variables are of particular importance for the analysis of the data from high 
energy heavy ion collisions. Furthermore, the calculation of the free energy by 
straightforward perturbation theory is infrared divergent already at three loops. 
This divergence, which is connected to the electric screening mass, can be re- 
summed, but at four loops new infrared divergences connected to the magnetic 
screening appear, which cannot be easily resummed. One may use dimensional 
reduction [H [21 El H], but the calculation of the free energy in this framework 
is a formidable task, which is not yet completely finished j5l El [7J. Therefore, 
at present a non perturbative theoretical method like lattice gauge theory is the 
only possibility to get quantitative results. This is particularly necessary below 
about twice the critical temperature, where perturbation theory even improved 
with dimensional reduction is not expected to work. 

The numerical demand of lattice QCD with fermions makes it difficult to make 
a continuum extrapolation with controlled systematical errors. For improved 
staggered fermions, however, there exist new data for the timelike lattice spacings 
N T = 4, 6, 8 [H [9]. In pure SU(3) gauge theory a continuum extrapolation has 
been performed a long time ago, but only up to T = 4T^|10j. In order to make 
contact with the different proposals for resummed perturbation theory, it would 
be very interesting to go to higher temperature. There are some recent results 
on this for N T = 4 with the Symanzik improved action 
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In this article we study a closely related theory, namely pure SU(3) gauge theory 
in 2 + 1 dimensions. This theory has many properties in common with 577(3) 
gauge theory in 3 + 1 dimensions. It has linear confinement and a deconfining 
phase transition at finite temperature. The infrared divergences in perturbation 
theory are stronger so that even the lowest non trivial order for the free energy 
is not calculable in straightforward perturbation theory [T2]. 

We determine the thermodynamical variables and the equation of state by lattice 
Monte Carlo computations on several lattices. In our calculations we use the 
so called integral method [TOj 13] . In order to change from the lattice coupling 
constant (3 to the temperature T, we use a scaling function derived from values 
for the string tension a at zero temperature [TU [03, [TBJ [T7j. Using the values 
for the critical couplings from Ref.[14j. we may also express the above thermo- 
dynamical variables as functions of T/T c . The method then gives directly the 
normalized trace of the energy momentum tensor, (e — 2p)/T 3 as a function of 
T/T c , where e is the energy density and p the pressure. Thermodynamical iden- 
tities in a sufficiently large and homogeneous system means that the normalized 
pressure, p/T 3 can be derived directly from the trace through integration over 
T. From these two quantities further thermodynamical variables, like the energy 
density e and the entropy density s can be easily calculated. We determine the 
thermodynamical variables in the range 0.6 < T/T c < 15, which in this the- 
ory means that the effective dimensionless coupling constant g 2 /T varies in the 
interval 3 > g 2 /T > 0.12. 

In Section 2 we present the problem and our adaptation of the integral method. 
In Section 3 we derive the /5-function and determine its parameters. In Section 4 
we present the results of our Monte Carlo calculation, and show that finite spatial 
lattice size effects are not present in the data. In Section 5 we discuss the results. 
Section 6, finally, is devoted to the conclusions. 

2 The method 

In this section we will first recapitulate some well known facts about thermody- 
namics and lattice gauge theory in 2 + 1 dimensions. This serves essentially to fix 
our definitions and notations. We then describe the method employed to extract 
the thermodynamical quantities from lattice data, following ref. P22 US] 
We start with the Euclidean Lagrangian in the continuum theory, 

C E (A^x)) = ^Ltr(F^(x)F^(x)). (1) 

Here, x = (xq,xi,Xz) is a three dimensional Euclidean vector. The dynami- 
cal variables are the gauge fields A^(sc), which are hermitean traceless matrices 
belonging to the algebra of the SU(3) group and 

F^(x) = d^A v (x) - d v A^(x) + i[A M (x), A v (x)\. (2) 
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Note that g 2 has the dimension of mass in three dimensions. 

The thermodynamics of the corresponding quantum theory is derived from 
the partition function, which is formally expressed as 

Z(T, V,g 2 ) = J VA^x)e~ti /Td *° I v d 2 xC E (A,(x)) ^ ^ 

where T is the temperature and V the spatial volume. 
From the partition function Z we get the free energy 

F(T,V,g 2 ) = -T\ogZ. (4) 

In the following we will use the pressure p, and the volume densities of the free 
energy, the internal energy and the entropy, and denote those by /, e and s 
respectively. We will assume that we have a large homogeneous system, in which 
case 

P(T, 9 2 ) = -f(T,g 2 ), (5) 
<T,g 2 ) = T 2 ^{p/n (6) 
s{T,g) = t^ = ^^. (7) 

It follows that in 2 + 1 dimensions we further have for the trace of the energy 
momentum tensor e — 2p, 

e — 2p d ( p 



J"3 



In fact we can use the left hand side as the basic quantity, from which the other 
thermodynamical quantities can be obtained apart from an integration constant. 

In 2 + 1 dimensions and for g 2 = (the free gauge theory), the temperature 
T is the only scale in the system, and the dimensionless quantities 

111 (9) 

are pure numbers. Thus, in this case the energy momentum tensor is traceless 
and 

6 = 2p, (10) 
. = f ■ (ID 

The pressure for eight free gluons in two spatial dimensions is easily calculated 
to be 

^ = 8^(3) = 1-5305.... (12) 



4 



In this article we consider the theory regularized on a finite lattice with lattice 
spacing a and with iV T points in the (inverse) temperature direction, defined as 
the direction and iVj points in the space directions 1,2. We denote the link 
variables starting from the site x in the positive directions by U^(x),fi = 0, 1,2 
and use the standard Wilson action: 

S w (U,(x)) = J>(f/ P ), (13) 
p 

S(U P ) = p(l-±ReTrUpY (14) 

where P denotes one of the 3N T x JV| plaquettes on the lattice and Up is the 
product of the [/-matrices around the plaquette. One should not confuse the 
dimensionless coupling constant (3 with the inverse temperature, which we will 
not use with that notation. Defining U^(x) in terms of A^x) by 

U^ x ) = e ii: +aii c* e ia ^ (x) , (15) 

one obtains the original Lagrangian from the classical limit of the lattice La- 
grangian, 



lim^F = ^TrF, v (x)F, v (x), (16) 



S(U P ) £ 

lim a(3 = 4- (17) 

Furthermore we define the temperature and the volume on the lattice by 

\ = aN T1 (18) 
V = (aN s ) 2 . (19) 

To obtain the thermodynamical quantities we need to calculate the free energy 
on the lattice, 

P = ~f = ^ogZ, (20) 
Z = / JJ dU^x)e- Sw(u ^ x)) . (21) 

The partition function cannot be calculated directly by Monte Carlo methods. 
Here we use the integral approach of Refs. [TQl E3] • Thus we first calculate the 
derivative of log Z with respect to the lattice coupling constant (3 at fixed iV r and 
Ns- It is given by 

-N T N 2 S (P s + 2P T ) . (22) 



dtegZ _ , x 

d/3 
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The negative sign may look strange, but we have not yet subtracted the pressure 
at zero temperature. Note that we define the plaquette average value (P) as 



{P) = (l-±ReTrU P \ (23) 

and Ps (resp. P T ) denotes a plaquette in the {1, 2} (resp. {0,1} or {0,2}) plane. 
Averages (. . .) are taken with respect to the Wilson weight corresponding to given 
lattice parameters. The averages containing either P T or P$ are always taken on 
the finite temperature lattices iV r < N$ and (Pq) will always denote average on 
the zero temperature lattice N T = N$. 

Eqs. (12"UI I2"2"j) defines p up to a constant with respect to /3, so that for some 
/3° to be chosen below we find 

rH 

a 3 p(/3, N T , N s ) = - dp (P s + 2P T ) + a 3 p((3°, N T , N s ). (24) 

We next subtract from this expression its value for iV T = Ns, which for Ns large 
enough constitutes its zero temperature limit, and thus consider the quantity 

a 3 (p(P,N T ,Ns)-p(P,Ns,N s ) 



fp0\* ( p o),3> - (Ps + 2P T )^ 



+ a 3 (p(p ,N T ,Ns)-p(P ,Ns,N s )). (25) 



The last term in this substracted expression for the pressure is still unknown, but 
it can be made negligible by choosing (3° in such a way that N T is sufficiently larger 
than the correlation length at that value of f3. To this extent then p(/3°, N T , Ns) is 
insensitive to the replacement of N T by Ns and the constant vanishes. Under this 
condition for (3°, our numerical estimate of the thermal (i.e. zero temperature 
subtracted) part of the pressure, which we denote by the same symbol p, is 

Pi(3, — = N " f d(3'AS((3', N T , N s ), (26) 

where 

AS((3, N T , N s ) = 3 (Pq)^ - (P s + 2P r ) /3 . (27) 

The above expression for the pressure is evaluated from simulations performed 
on N T x Ng lattices, with the aspect ratio Ns/N T > 8. In section 4 we will show 
that for this minimum aspect ratio we are sufficiently close to the thermodynamic 
limit, i.e. that the finite size corrections are smaller than our statistical errors. 
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In order to obtain the pressure from (126]) as a function of T/T c , where T c is 
the critical temperature of the continuum theory, we still need to relate T/T c to 

ft 

(3 = (3(T/T C ), (28) 

which will be discussed in Section 3. Finally we choose (3° = /3(T/T C = 0.6), after 
checking from the measurements that iV 3 times the integrand in (126!) is negligibly 
small at this temperature. According to (18|26II28I) . the normalized trace of the 
energy momentum tensor is finally given by 

e - 2 P T _d_ (P_\ 
T 3 dT VT 3 / 

= iV 3 AS^(|),iV T ,iV^ T^. (29) 

From Eqs. (171), (|26|) and fl29j) one obtains the expressions for the energy and 
entropy densities e and s. 



3 Physical Scales for the Temperature 

We need to express the thermodynamical variables, provided by (|26|) and (l29j) 
from lattice simulations, as functions of the lattice temperature T defined in (ITS]) . 
This temperature has to be measured at some scale related to the continuum 
theory. The continuum coupling g 2 is such a scale, but it is interesting to consider 
other quantities of the same dimension, and thus proportional to g 2 , which have 
a direct physical interpretation. Here we will use a/cx, the square root of the 
zero temperature string tension, and T c , the critical temperature, both in the 
continuum theory. 

In this section, we construct the corresponding scaling functions which, given 
Ns and N T , represent the change of variable from (3 to T/y/a, T/T c or T/g 2 . We 
start with the former, namely 

P = f3 a {T/^) and its inverse -= = (30) 

Passing from one scale to the other is just a multiplicative numerical factor, 
T c j y/a, which must be also computed. 

In the literature, there exist mainly two sets of data for the quantity 

F a {(3) = a^, (31) 

obtained from lattice simulations in various domains of Ns and (3. A large set 
(3 G [8 : 50], referred to as A in what follows, has been covered in [151 EE], 
while recent simulations [17] . provide high accuracy data at three (3 values only, 
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24 
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0.20712(24) 




19 
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0.19447(23) 
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21 




0.173948(75) 
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22 


32 


0.16555(24) 


na 


24 


32 


0.15112(23) 
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28 


32 


0.1275(2) 


[16j 


34 


40 


0.10379(26) 


[16j 


40 




0.087046(75) 


m 


50 


48 


0.07021(5) 


[15j 



Table 1: Data for the string tension. Each line gives the values of (3, Ns, and 
F a {(3) — a an d the corresponding reference. 

14.717, 21 and 40 (set B). Both sets are included in table [U How we use them 
is explained below. 

For any value of (3 occuring in Table [fl using Eqs. (I31II18P gives directly 

T 1 

7^ ~ N T F a {py 

However, we need in practice to evaluate T/y/a{f3) for a denser and larger set 
{(3 G [8 : 400]) of (3 values. This is required in particular to obtain the pressure 
from the integral (T26"]) . in a range of temperatures T/T c from about 0.6 to 15. 
This we achieve by making fits to the data of Table [Q We use the form 

F -e> = | + ? + ?- < 33 > 

This choice is done for the sole purpose of representing properly the data, with 
no theoretical prejudice about the actual behaviour of F a around l//3=0. An 
investigation shows that there exists a systematic although small discrepancy 
between the sets A and B, and we found it adequate to fit them separately 
rather than altogether. Replacing the symbol q of (13"3"j) by aj and 6j respectively 
for A and B, the two fits give 

ai = 3.37(l) a 2 = 3.90(25) a 3 = 50.1(1.8) (34) 
6i = 3.34(1) b 2 = 4.68(50) b 3 = 39.8(5.0). (35) 
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The relative difference of 1% between a\ and b\ gives the order of magnitude of 
the systematic effect due to choosing A or B at large (3. At low (3, one might 
worry that the large cubic term found in both cases casts a doubt on the validity 
of the extrapolations. We note however that table [T] gives F a = 0.568(2) at 
(3 = 8.156 [15], less than 3% away from the average values 0.564 and 0.553 one 
gets respectively (with larger errors) from A and B. Furthermore, as we shall see 
in the next sections, such low values of (3 are used only for N T = 2, which plays 
a very marginal role in our main physical results. With that in mind and for 
definiteness, we will use the parameters flHBT) corresponding to the high precision 
data [17J. The numerical value of T j \fo given /3, N T is thus obtained by inserting 
(|33|) in (|32|) with q = fej, bi as in (135]) . The errors on yfa in the input data imply 
an error on this change of variable, to be estimated (see below). 

Our last task before that is to define T using the scale T c , the critical temper- 
ature in the continuum, that is in the limit N T and N$ — > oo. Since we already 
know T/y/a we thus compute T c j 1 \fo. From fl32|) . we write 

yfd Nt-*oo N T F C r[ft c (N T )Y 

where j3 c {N T ) is the critical lattice coupling determined by lattice simulations at 
given N T . In this expression, the N$ = oo limit should have been taken first. 
As will be shown in section 4 on the basis of our simulations, the value of N$ 
(8xN T ) for each N T is large enough to eliminate any finite N$ effect. In order to 
perform the limit fl36l) . we use the set of (3 C obtained in [15] for N T —2, 4 and 6, 
that is respectively /3 C =8. 155(15), 14.74(5), and 21.34(4)(11). More recent data 
on (3 C {N T ) at N T = 2,3,4,5 can be found in p2] and they do not modify our 
estimate, which is: 

T 

— = = 1.00(4). (37) 

The ratio to this constant of T j \fa is T/T c , which will be used as the temperature 
variable in the rest of the paper. 

The error in (1371) , as well as any other error quoted in the present and subsequent 
sections have been estimated using the so-called bootstrap method (see [T9l 120]). 
To be concrete, we describe it in detail for this very typical case. 

The set of initial data, noted Xj for short, consists of values and errors for 
a \fo and (3 C at given values of external variables (3 or N T . All Xj are assumed to 
be independent and gaussians. We performed the following steps 

Step 1: From the last assumption, generate another equivalent set of input 
variables jjj. 

Step 2: Fit the coefficients q of (1331 to the yj values. (Note that in the case 
where yj represents the data of [17] . the 3 coefficients aj are exactly determined 
by the 3 yj. Standard \ 2 methods would not give any estimate of their errors). 

Step 3: Compute F a {(3) for (3 = (3 C (N T ) 
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(36) 



Step 4'- Fit a quadratic form in 1/N T to the 3 values obtained for the inverse 
of N T F a (/3 C (N T )) . From Eq. (13Tj) . its value at 1/N T =0 is an estimate of T c /y/a. 

Step 5: Repeat everything from Step 1 many times (100), each time with a 
new set of jjj. 

Our result fl37j) gives the value of T c j ^fa for = Xi, with an error equal to 
the standard deviation measured within the set of estimates obtained at Step 5. 



A different way to determine a (3 function of the form (130]) is to use directly 
the knowledge of (3 C available for various iV T values. We expand (3 as a Laurent 
series in ag 2 , writing 

6 

P = —= + b + cag 2 + (38) 
ag 2 

The first term corresponds to the classical limit of the lattice action. We include 
as many terms as there are known (3 C values, and consider the truncated series 
obtained as representing the first terms of an asymptotic series. We then replace 
the lattice spacing a by 1/ (N T T) and introduce the scale T c to rewrite (1381 as 

P 6T C T b g 2 T c 1 

This expansion now constitutes an expansion in 1/N T , whose 1/N T = limit 
gives an estimate of T c in units of g 2 in the continuum limit. 

From Eq. f[3"9~j) and the values of P c quoted in [T5] for N T = 2, 4, 6, namely 
8.149(3), 14.74(5), 21.34(11) respectively, one finds 

T 

-4 = 0.55(2); 6=1.5(6); c = 0.06(50). (40) 
9 

We observe that the 1/N 2 term is compatible with zero in the range of interest, 
which means that the P function in terms of T/g 2 is essentially linear. Although 
the P functions determined via the knowledge of the string tension and via Eqs. 
fl3"9"| |4"U|) need not be the same at finite iV T , we find them very close one to the 
other, so that they can be used equivalently in the present work. In practice, 
an accurate representation of our findings for the /^-function, valid in the whole 
range explored in P and N T is 

^ = 3.3^ + H. (41) 

This formula is useful in order to jump easily from functions of P to functions 
of T/T c . Only when precise estimates of errors in final results is needed is it 
necessary to go through the process described above and apply the bootstrap 
method. 
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4 Results of the Simulations 



The results of our simulations consist of high precision plaquette data, that is 
their average values and statistical errors, for zero and finite temperature lattices. 
These data are shown in the appendix Here we present and discuss the quantity 
N% AS = N 3 (3 (Pq) — (2P T + Ps)) which was shown in section 2 to constitute 
the basis for all computations of the thermodynamical variables. 

An important source of possible systematic errors for the quantities of physical 
interest resides in the finite size effects (FSE) at (3 and N T fixed due to Ns being 
finite. A reason for this is that although (TrUp) / 3 is of order one, very strong 
cancellations occur in AS. According to the definitions in section 2, (P) is 
expected to go to as (3 — > oo. Although its precise behaviour around = 
is unknown, we assume that it can be numerically represented by a (probably 
asymptotic) series in We find that the linear term is absent and that N% AS 
is of the form 

N 3 T AS « A + (/r 3 ), (42) 

at large /3, with A a function of N T and Ns- This behaviour was verified to a 
high degree of precision from our data. 

We, of course, have to make sure that the finite size effects occuring in the 
distinct measurements of Pq,P t , Ps are so small that the corresponding effects in 
N% AS lead to acceptable errors, smaller than the statistical ones. These finite 
size effects should be related to the largest gauge invariant correlation length of 
the system, a priori different for P , and for Ps and P T . In the former case we 
are in the confined phase at T = 0, where the largest correlation length is given 
by the inverse of the lowest glueball mass, whose measurement was reported in 
|16j . In the latter case, it is given by the screening length, measured in [21 J . 

Here we report the results of a systematic study of the FSE observed in 
P , P T , Ps, from simulations performed, at N T =4 fixed, for a large set [15:150] of 
f3 values and N s e [6 : 48]. 

Each (P) at a given (3 was analyzed as a function of Ns, and parametrized as 

{P) PtNB = A(J3) - B{(3) eM~KP)N s ). (43) 

For a lattice spacing a, this corresponds to the assumption that an effective 
correlation length 

m = a//i(/3) (44) 

leads to FSE which decay exponentially in the ratio L/£, where L = aNs is the 
spatial extent of the lattice. By fitting Eq. (jl"51) to the data, we determined A(/3), 
the plaquette value on a spatially infinite lattice, the weight B((3) of the exponen- 
tial component, and the effective mass fi{(3) in lattice units. Depending on the 
type of plaquette considered, the quantities A,B,fi are given the corresponding 
subscript or r or S. 
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Figure 1: The average spatial plaquette (Ps(/3,Ns)) approaches its limit 
(Ps(P, oo)) from below, exponentially in N$. For clarity of the figure, the quan- 
tity A s {f3) - (P S (/3,N S )) has been multiplied by (3 2 . 

The general features of our findings on these FSE are the followings. 

1. The parametrization ( )43l) is always adequate. As an illustration, Fig. (TjQ) 
shows how the data for (Ps(fl, N s )) approach their N s = oo limit, together 
with fits to Eq. (j43l). 

2. Close to P c , the B coefficients are of opposite signs in P T and P$, and a 
partial cancellation of the FSE occurs. They are of the same positive sign 
in P Q and Ps, which again lowers the overall FSE. 

3. The mass parameter fi T , found to be around .45 for (3/N T large is poorly 
determined. On the contrary, /is is well determined, and it lies in between 
the two screening masses Ms and Mp found in [21], as illustrated in Fig. §2§. 
In the range (3 & [15 : 34] for which data on glueball masses exist [IE] , we 
find that /i is lower than those, but only by about 20%. In the whole range 
/3 E [15 : 150] we have approximate scaling, i.e. (3 Ho{(3) oc /i (/5)/a is nearly 
constant. 

4. Finally the comparison with the statistical errors of the effects due to Ns 
finite shows that the latter can be neglected in our final measurements, all 
performed with the aspect ratio £ = Ns/N T = 8. In the present exploration, 
we may compare the cases £ = 4 and £ = 8 corresponding respectively to 
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Figure 2: 

Comparison of the inverse correlation length [is appearing in the FSE for P$ with 
the screening masses Ms and Mp determined in Ref. |21j. 

Ng=16 and 32. Then for each (P), the AT 5 — dependent term in (|43p is 
reduced by a factor exp(16 /i) if £ = 8, that is at least by two orders of 
magnitude in typical cases. Moreover, as mentioned in point 2 above, the 
contributions of these terms to AS from different plaquettes tend to 
cancel each other. On the contrary the statistical errors are much larger in 
N% AS than in (P) in relative value. 

The results for N% AS, for those lattice sizes of particular interest in the final 
computation of the thermodynamical variables (iV T = 4, 6, 8) are presented in 
figure [3l Their numerical values are given in the tables of the appendix, where 
the case N T = 2 is also included. 



5 Thermodynamics 

From the results presented, just above for N% AS as a function of the lattice 
coupling (3 and in section 3 for the scaling function relating (3 to the temperature, 
the three thermodynamical variables e, p and s can be explicitly determined from 
the formulae of section 2. 
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Figure 3: The average plaquette combination N^AS for N T = 4, 6, 8. Their 
values and errors (in most cases smaller than the point sizes) are tabulated in 
the appendix. On each plot the horizontal scale is (5 / N T for easier comparison. 
The vertical scale is the same on all plots. 
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For definitness, we choose to measure the temperature in units of the critical 
temperature. The temperature variable and the scaling function are then denoted 



t = ^, (45) 
(5 = P(t). (46) 

On a given N T x iVj lattice, the dimensionless energy and pressure densities 
e/T 3 and p/T 3 follow from Eqs. (I26| 129]) rewritten as 

( ~ 2P it) = N 3 As(p(t))A (47) 



T 3 r V w / dt 



p(t) 
T 3 



iV 3 / dp' AS 03'). (4* 



In this integral, we choose to — 0.6, a temperature below which the integrand is 
always negligible. From these equations one directly obtains p/T 3 and e/T 3 ; the 
entropy density s/T 2 follows from fl7|). 

In order to use Eqs. (I47p in practice, and in particular to perform the in- 
tegral defining the pressure conveniently, we have interpolated N 3 AS(/3(t)) via 
a smooth function f(t). In the absence of any theoretical ansatz for the shape 
observed in figure [HJ we choose the simplest parametrization adapted to the de- 
scription of a sharp rise in the vicinity of the critical point t—1, and the t~ 2 fall 
off at large t expected from Eqs. (l42|4ip . The following function constitutes an 
accurate representation of N 3 AS((3(t)) within errors 

f(t) = ait i + a3e (49) 

The adequacy of this function for our purposes is illustrated by the curves 
drawn through the data points in figure HI The horizontal error bars represent 
the errors in the change of variable f3 — > t = T/T c described in section 3. They 
are of the order of 3 percents of t. 

One may worry that choosing a specific function like (j4"9~i) introduces a sys- 
tematic error into subsequent calculations. A different choice was tried during 
this investigation. To the extent that it provided a comparable agreement with 
the data to which it was fitted, it leads, for example for the pressure, to differ- 
ences much smaller than the statistical errors. However one must keep in mind 
that although f(t) constitutes a good interpolation of the data within errors, the 
systematic errors attached to its use above t ~ 15 are not controlled. 

Our final results for the pressure at iVy = 2, 4, 6, 8 with Ng = 8 iV T are shown 
in figure [5j The errors on the scale t, which we showed in figure HJ do not affect the 
pressure, which varies slowly where they become important (at large t). For the 
largest lattice considered (iV T = 8), figure [6] compares the three thermodynamical 



15 




Figure 4: Trace of the energy momentum tensor. 



variables, after multiplication by coefficients such that they would be the same in 
the free theory (see Eq. [TO]) . The continuum Stefan-Boltzmann value for a free 
gas (horizontal dotted line) is also shown. 

All the errors on the quantities plotted in figures HJ and E] have been esti- 
mated using the bootstrap method, similarly to the case described in section 3 
for the calculation of T / ^/a. 



Let us now comment these results, whose extensions towards the continuum 
limit 1/N T = and to larger temperatures are highly desirable. A look at the 
results suggests that the continuum limit is nearly reached at N T = 8, the more 
so t is large. We observe that the values of (e — 2p)/T 3 , at least above t ~ 5 
hardly distinguish between N T = 6 and 8, and that the difference between the 
corresponding pressures does not exceed ~ 2a. More quantitative statements 
however are hard to produce. Attempts at fitting quadratic functions in 1/N T or 
1/iV 2 to the pressure at fixed t remained inconclusive, with no strong constraint 
on the continuum limit. 

We then turn to the question of the high temperature behaviour. We already 
noticed (02} that the tail at large (3 of iV 3 AS was compatible with l//? 2 , which 
implies via (j4"Ti) 

e - 2p b 



T 3 t 



t > ti. (50) 
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Figure 6: The pressure, and the energy and entropy densities for the largest 
lattice size (8x64 2 ). 
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train 


b 


P/T (train) 


p/T 3 (oo) 


4.0 


1.81(7) 


1.13(2) 


1.58 


6.0 


1.7(1) 


1.27(2) 


1.56 


8.0 


1.6(2) 


1.34(2) 


1.53 



Table 2: 



for some large enough t\. Since this appears to be true for any iV T , it is reasonable 
to expect that it is also valid in the continuum. We thus formulate the following 
ansatz: For t larger than some t m i n , the continuum limit of (e — 2p)/T 3 is close 
to its value measured at N T = 8 and it can be parametrized there by b/t. We 
then obtain from Eqs. (j4"Tl |4T|) 



> to 

b 

= a — -. 

f 



(52) 



a = if^(train) + ~ • (53) 

t m in 

For N T = 8, the data between t ~ 4 and t — 15 (the largest value explored) is 
compatible with a behaviour of the type fl50|) . Suppose now that this remains 
true above £=15. Then from f)52ll53p . we have 

j$(t = oo) = ^(U) + —. (54) 



Given t min , we get b from (e - 2p)/T 3 (t m in), and p/T 3 (t m in) is known from jSTJ), 
so that a consistency check is that p/T 3 (t = oo) is independent of t m j n . In table 
[21 we report the values found for b and p/T 3 (t = oo) for t min = 4,6 and 8. We did 
not try to estimate any error on p/T 3 (t = oo), certainly affected by uncontrolled 
systematics. Its value is stable with respect to t min , and although it is obtained 
for N T finite, it is found close to the continuum Stefan-Boltzmann value for a free 
gas, 1.53. 



6 Conclusions 

In this article we have determined the thermodynamical quantities for SU(3) 
gauge theory in two spatial dimensions. This determination has been performed 
through lattice simulations, measuring the thermodynamical variables as func- 
tions of the lattice coupling j3. We then constructed the scaling function connect- 
ing j3 to the temperature expressed in units of various physical quantities, the 
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square root \fa of the zero temperature string tension, the coupling constant g 2 
of the continuum theory and the critical temperature T c . In particular, we deter- 
mined the ratios T c j ' \fa and T c j ' g 2 to be respectively 1.00(4) and 0.55(2). This 
construction has been discussed in detail, using published data on the string ten- 
sion and on the critical lattice coupling f3 c .We finally kept T c as the temperature 
scale and presented the thermodynamical variables as functions of T/T c . 

We have determined the finite volume N$ effects at fixed bare coupling con- 
stant (3 by a careful analysis and shown that they are negligible compared to 
the statistical errors. We have further shown that those finite size effects fall off 
exponentially with N$, and that the coefficient in the exponent multiplying N$ 
is consistent within about 20% with the longest correlation lengths expected on 
the various lattices, namely the inverse of the lowest glue ball mass on the zero 
temperature lattices, and the lowest screening length on the finite temperature 
lattices. This has lead us to choose the aspect ratio N$/N T = 8, for which the 
finite size effects are much smaller then the statistical errors. 

We investigated the finite lattice spacing effects by simulating lattices with 
extension in the temperature direction N T = 2, 4, 6, 8. We found that N T = 8 is 
large enough to be near to the continuum limit. We have not presented attempts 
to extrapolate the thermodynamical variables to the continuum, because addi- 
tional systematic errors appear. Instead we have analyzed the thermodynamical 
variables for N T = 8. We find that for large enough temperatures, well above 
the phase transition region, the pressure can be represented by a very simple 
function, namely p = T 3 (a — bT c /T). Using the thermodynamical identities for 
a large homogeneous system we obtain the corresponding formulae for the other 
thermodynamical quantities. 

Improving our findings in view of a quantitative determination of the con- 
tinuum limit requires a more precise analysis of the N T dependence of the basic 
quantity N 3 AS. Larger statistics and/or larger lattices are possible ways in this 
direction. 
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A Data 

The data were produced from a standard SU (3) MC program. We used a mix- 
ture of heatbath and overrelaxation moves, with typically one heatbath sweep 
after 2 to 4 overrelaxation sweeps. All the simulations were done on the COCOS 
opteron cluster in the Institute of Physics, Jagiellonian University. This ap- 
pendix presents the data for the plaquette averages (P x ) — 1 — g Tr (Up x ). The 
subscripts S and r denote spacelike and timelike plaquettes respectively, on finite 
temperature lattices (N T < Ng). The subscript zero denotes plaquettes on zero 
temperature lattices (N T = Ng). For each of the N T x Ng or N s lattices, as 
functions of j3 we list the values of (P T ), (Ps), (Po), and (N^AS) from which 
the thermodynamical variables are finally computed. The errors were estimated 
using standard methods. We took care of strong correlations between P T and 
Pg, by first measuring the statistical errors for the combination 2P T + Pg, and 
adding them quadratically to those of the independent quantity 3Po to finally get 
the errors on (AS) = (3P - (2P T + Pg)) 
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Table 3: Plaquette data for lattices of size 2xl6 2 and 16 3 . 
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(Pr) 
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Table 4: Plaquette data for lattices of size 4x32 2 and 32 3 . 
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U.lozz000(zo J 


n 1 90999/2/90^ 

U.lozooo0(ozJ 


n 1 909901 /i/;^ 
U.loZooUl(10 J 


n 0901 /i o\ 
U.UoUl(lo ) 


oo nnn 
zz.UUU 


U.lzOyoUo(zO ) 


n 1 okoc;90/91 "\ 
U. lz0U0oU(ol ) 


n 1 ofin97fi/i r\ 
U.lzOUo /0(10 ) 


n n/i9c;/i7 N \ 
U.U4oO(l / ) 


no nnn 
zo.UUU 


n 1 om noi /1 a\ 
U.lzUlyzl(10J 


n 1 on99nn/i c"\ 
U.lzUooyU(loJ 


n 1 0091 009/C9\ 

U.lzUolyzo(oo J 


n n/inoo/7n"\ 
U.U4yzz( /y ) 


nc nnn 
zO.UUU 


U.1U0 (o (oyZO ) 


U. lUOyol l(zo ) 


n 1 nt^m ^1/^/7/1^ 
U.lU0yl010( /4J 


n n/i /i 9fi /7o"\ 
U.U44oo( / y J 


on nnn 

zy.uuu 


U.Uy44olO(l I ) 


n nn /i c o /oo\ 
U.Uy400o2(zU J 


n nn /i ^noo^ /£; tr \ 
U.Uy40U2o0(0o ) 


n n/t f\A /i n\ 
U.U4U4(1U J 


33.000 


0.0827075(15) 


0.0828960(18) 


0.08281948(71) 


0.03179(92) 


37.000 


0.0735500(13) 


0.0737324(16) 


0.07365035(80) 


0.02530(86) 


41.000 


0.06622256(75) 


0.06639958(91) 


0.06631278(43) 


0.02032(47) 


48.000 


0.0563954(11) 


0.0565549(13) 


0.05647074(47) 


0.01462(63) 


60.000 


0.04495852(65) 


0.04510118(80) 


0.04502038(34) 


0.00925(39) 


80.000 


0.03360653(36) 


0.03372051(45) 


0.03365228(25) 


0.00519(24) 


100.000 


0.02683241(36) 


0.02692816(44) 


0.02686872(16) 


0.00304(21) 


120.000 


0.02233119(22) 


0.02241268(28) 


0.02236182(17) 


0.00223(15) 


150.000 


0.01784145(19) 


0.01790881(24) 


0.01786618(13) 


0.00152(13) 


200.000 


0.01336444(13) 


0.01341620(17) 


0.01338281(10) 


0.000763(91) 


240.000 


0.01112982(12) 


0.01117361(15) 


0.011145438(75) 


0.000661(77) 


300.000 


0.00889827(10) 


0.00893393(13) 


0.008910796(50) 


0.000391(62) 



Table 5: Plaquette data for lattices of size 6x48 2 and 48 3 . 
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p 


(Pr) 


(Ps) 


(Po) 


(N* AS) 


10.000 


0.2975884(59) 


0.2976005(72) 


0.2975867(40) 


0.0092(60) 


20.000 


0.13928938(92) 


0.1392889(12) 


0.13929140(79) 


0.0028(25) 


23.000 


0.12031707(79) 


0.1203197(10) 


0.12031754(73) 


-0.0019(21) 


25.000 


0.11031138(86) 


0.1103157(11) 


0.11031518(72) 


0.0054(19) 


27.000 


0.10184029(88) 


0.1018568(11) 


0.10185383(82) 


0.01591(71) 


29.000 


0.09456979(81) 


0.0946071(11) 


0.09460236(76) 


0.03296(65) 


30.000 


0.09131473(79) 


0.0913588(10) 


0.09135313(72) 


0.03632(72) 


31.000 


0.08828064(88) 


0.0883287(11) 


0.08831814(80) 


0.03708(61) 


33.000 


0.08277950(54) 


0.08283196(70) 


0.08281902(50) 


0.03451(62) 


35.000 


0.07792616(50) 


0.07798094(66) 


0.07796591(47) 


0.03257(58) 


37.000 


0.07361228(63) 


0.07366982(82) 


0.07365076(58) 


0.03001(55) 


40.000 


0.06796965(58) 


0.06802736(76) 


0.06800643(53) 


0.02589(46) 


50.000 


0.05414363(41) 


0.05419933(54) 


0.05417403(37) 


0.01708(44) 


60.000 


0.04499613(38) 


0.04504689(50) 


0.04502046(35) 


0.01166(32) 


80.000 


0.03363355(33) 


0.03367621(43) 


0.03365276(30) 


0.00680(32) 


100.000 


0.02685397(17) 


0.02689045(22) 


0.02686885(15) 


0.00436(23) 


120.000 


0.02234902(13) 


0.02238078(17) 


0.02236201(12) 


0.00290(18) 


140.000 


0.01913902(16) 


0.01916675(21) 


0.01914968(15) 


0.00216(16) 


160.000 


0.01673517(14) 


0.01676024(19) 


0.01674472(13) 


0.00172(14) 


180.000 


0.01486781(12) 


0.01489035(16) 


0.01487629(11) 


0.00139(12) 


200.000 


0.01337535(11) 


0.01339589(15) 


0.01338306(10) 


0.00102(11) 


250.000 


0.010692216(90) 


0.01070908(12) 


0.010698237(84) 


0.000777(90) 


300.000 


0.008905754(75) 


0.008920037(99) 


0.008910804(69) 


0.000474(76) 


350.000 


0.007630792(65) 


0.007643064(84) 


0.007635106(60) 


0.000445(65) 


400.000 


0.006675004(80) 


0.00668602(10) 


0.006678986(93) 


0.000243(59) 



Table 6: Plaquette data for lattices of size 8x64 2 and 64 3 . 
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